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Abstract 

We consider two quantum cosmological models with a massive 
scalar field: an ordinary Friedmann universe and a universe contain- 
ing primordial black holes. For both models we discuss the complex 
solutions to the Euclidean Einstein equations. Using the probability 
measure obtained from the Hartle-Hawking no-boundary proposal, we 
find that the only unsuppressed black holes start at the Planck size 
but can grow with the horizon scale during the roll down of the scalar 
field to the minimum. 
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1 Introduction 



In this paper we ask how likely it is for the universe to have contained 
primordial black holes. We investigate universes which undergo a period of 
inflation in their earliest stage, driven by a scalar field <fi with a potential V(<f>) 
with a minimum V(0) = 0. The results do not depend qualitatively on the 
exact form of the potential, so for simplicity we consider a massive minimally 
coupled scalar V((f>) = |m 2 2 . The scalar field starts out at a large initial 
value (fro and acts as a cosmological constant for some time until it reaches the 
minimum of its potential and inflation ends. We consider two different types 
of spacetimes: in the first, the spacelike sections are simply 3-spheres and 
no black holes are present; in the second, they have the topology S 1 x S 2 , 
which is the topology of the spatial section of the Schwarzschild-de Sitter 
solution. Thus these spaces can be interpreted as inflationary universes with 
a pair of black holes. In the inflationary period, the first type will be simi- 
lar to a de Sitter universe, the second to a Nariai universe [H. To find the 
likelihood for primordial black holes, we assign probabilities to both types of 
spacetimes using the Hartle-Hawking no-boundary proposal (NBP) [[J. This 
is the only proposal for the boundary conditions of the universe that seems 
to give a well-defined answer in this situation. It is not clear how to apply 
the so-called 'tunneling proposal' in the S l x S 2 case. If one takes the action 
to appear with the opposite sign as is done in the S 3 case, one would reach 
the conclusion that a universe with a pair of black holes was more likely than 
a universe without, and that the probability would increase with the size of 
the black holes. This is clearly absurd. 

The NBP framework is summarized in Section [2|. In Sections |3| and ^ 
we review its implementation for cases with a fixed cosmological constant. 
In Section [| we introduce a massive scalar field and discuss the solutions 
of the Euclidean Einstein equations for the S 3 case. They will be slightly 
complex due to the time dependence of the effective cosmological constant 
(m(j))~ 2 . We obtain the Euclidean action for those solutions. In Section |5] 
we go through a similar procedure for the S 1 x S 2 case. We find that the 
black hole grows during the inflationary period, a noteworthy difference to 
the Nariai case with a fixed cosmological constant. In Section [7] we use the 
action to estimate the relative probability of the two types of universes. We 
find that black holes are suppressed for all but very large initial values of 0o- 
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2 The Wave Function of the Universe 



The Hart le- Hawking no-boundary proposal states that the wave function of 
the universe is given by 

tfofoj, $9m] = J D( 9lxu , $) exp [-I{g^ $)] , (2.1) 

where (/i^, $aAf) ar e the 3-metric and matter field on a spacelike boundary 
dM and the path integral is taken over all compact Euclidean four geometries 
g^y that have dM as their only boundary and matter field configurations $ 
that are regular on them; I(g IM/ , 3>) is their action. 
The gravitational part of the action is given by 

^ = -77T- / d A xg l /\R- 2A) - J- / d 3 */* 1 ^ (2.2) 

lD7T JM 87T JdM 

where R is the Ricci-scalar, A is the cosmo logical constant, and K is the 
trace of Kij, the second fundamental form of the boundary dM in the metric 
g. For the origin of the boundary term, see, e. g., ref. |§. 

In the standard 3+1 decomposition M, the metric is written as 

ds 2 = N 2 dr 2 + h^dx 1 + N i dr)(dx j + N j dr). (2.3) 

Assuming that the NBP is satisfied at r = 0, the Euclidean action then takes 
the form 

I E = — — I™' Ndr [ d 3 x h^t-KalC 3 + K 2 + 3 R- 2 A) 
16n J T =o J 

+ — [ d 3 xh l ' 2 K. (2.4) 
8n Jt=o 

Here 3 R is the scalar curvature of the surface, and tensor operations are 
carried out with respect to the surface metric hij. In the first term the 
boundary terms are implicitly subtracted out at r = and r = tq M . But it 
is an essential prescription of the NBP that there is no boundary at r = 0. 
So the second term explicitly adds the contribution from r = back in. 
It vanishes for universes with spacelike sections of topology S* 3 , but can be 
non-zero for the topology S 1 x S 2 . 
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There are unresolved questions on how to choose the integration con- 
tour and make the integral converge ||, but we shall not discuss them here. 
Instead, we will use the semiclassical approximation 



-In 



(2.5) 



where the sum is over the saddlepoints of the path integral, i. e. the solutions 
of the Euclidean Einstein equations. In this paper, we neglect the prefactors 
A n and take only one saddlepoint into account for a given argument of the 
wave function. So the probability measure will be 



dM 



_ 2/ Rc 



(2.6) 



where 7 Re is the real part of the Euclidean saddlepoint action. 

By considering only spaces of high symmetry (homogeneous S 3 or S 1 x S 2 
spacelike sections) we restrict the degrees of freedom in the metric to a finite 
number q a . The Euclidean action for such a minisuperspace model with 
bosonic matter will typically have the form 



Ndr 



1 dftrf 



(2.7) 



The saddlepoints will in general be complex solutions q a (r) in the r-plane. 
In the semiclassical approximation the following relations for the real and 
imaginary parts of the saddlepoint actions hold: 



I( V / Rc ) 2 + I(v/ Im ) 2 + t/(g a )=0 



V/ Re • V/ 



im 



0. 



(2.8) 
(2.9) 

where the gradient and the dot product are both with respect to f a @. There- 
fore J Im will be a solution of the Lorentzian Hamilton- Jacobi equation in 
regions of minisuperspace where \1/ has the property that 



(w 



^ 2 <{:: fv/ Im ^ 2 



(2.10) 



This allows us to reintroduce a concept of Lorentzian time in such regions: We 
find the integral curves of V/ Im in minisuperspace and define the Lorentzian 
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time t as the parameter naturally associated with them. Reversely, if we 
demand that the NBP should predict classical Lorentzian universes at suffi- 
ciently late Lorentzian time, condition ( |2.10p must be satisfied. This means 
that there must be saddlepoint solutions for which the path in the r-plane 
can be deformed such that it is eventually almost parallel to the imaginary 
r axis and that all the q a should be virtually real at late Lorentzian times. 
In summary, the following conditions must be met: 

i. The NBP must be satisfied at r = 0. 

ii. At the endpoint tqm of the path, the q a must take on the real values 
qg M of the arguments of the wavef unction: 

q a (r dM ) = q a 9M - (2.11) 



iii. The q a must remain nearly real in the Lorentzian vicinity of the end- 
point: 

' dq° 
dr 



Re 



0. 



(2.12) 



3 The de Sitter Spacetime 

In this and the next section we review vacuum solutions of the Euclidean Ein- 
stein equations with a cosmological constant A. First we look for a solution 
with spacelike sections S 3 . Therefore we choose the metric ansatz 

ds 2 = N(r) 2 dr 2 + a(T) 2 dnl (3.1) 

The Euclidean action is 

3ir f , / a 2 A 9 \ , , 

I=--jNdra[^- 2 +1-3« 2 J, (3-2) 

A dot (') denotes differentiation with respect to r. We define 
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Variation of a and N yields the equation of motion 



a ~ T 2 



+ H = (3.4) 
a 



and the Hamiltonian constraint 

a 2 1 



-2 - -o + # = (3.5) 
er er 



in the gauge A/" = 1. A solution of equations ( |3.4| ) and ( |3.5| ) is given by 

a(r) = H- 1 sin Ft. (3.6) 

It is called the de Sitter spacetime. The NBP is satisfied at r = 0, where a = 
and 4^ = 1. If we choose a path along the r Re -axis to r = ^ , the solution 
will describe half of the Euclidean de Sitter instanton S* 4 . Choosing the path 
to continue parallel to the r Im -axis, a(r) remains real and the conditions i to 
Hi of the previous section will be fulfilled: 



-Rc— 7T 

2H 



H- 1 cosh #r Im (3.7) 



This describes half of an ordinary Lorentzian de Sitter universe. 

So with the above choice of path, equation (|3.6|) corresponds to half of a 
real Euclidean 4-sphere joined to a real Lorentzian hyperboloid of topology 
R 1 x S 3 . It can be matched to any clq M > by choosing the endpoint 
appropriately, and for ciqm > H~ l the wavefunction oscillates and a classical 
Lorentzian universe is predicted. 

The real part of the action for this saddlepoint is 

Cs itter = | J* dr^a (H 2 a 2 - l) = ~ (3.8) 
The Lorentzian segment of the path only contributes to 7 Im . 



4 The Nariai Spacetime 

We still consider vacuum solutions of the Euclidean Einstein equations with 
a cosmo logical constant, but we now look for solutions with spacelike sections 
S 1 x S 2 . The corresponding ansatz is the Kantowski-Sachs metric 

ds 2 = N(r) 2 dT 2 + a(r) 2 dx 2 + b(r) 2 dtt 2 2 . (4.1) 
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The Euclidean action is 



I = —7i / Ndra 



'ii 

N 2 



b ab 
''aN 2 



1-Ab 2 



7Y 



—ab — labb 



T=0 



(4.2) 



where the second term is the surface term of equation (|2.4ft . We define 

H = VX. (4.3) 

Variation of a, b and iV gives the equations of motion and the Hamiltonian 
constraint: 



ab 
ab 



b 
b 

2— - 

ab 



ab a tt9 
- — + - + H 2 
ab a 





0. 



A solution is given by 



(4.4) 
(4.5) 
(4.6) 

(4.7) 



a(r) = H 1 sin Ht, b(r) = H 1 = const. 

It is called the Nariai spacetime. The NBP is satisfied at r = 0, where 

a = 0, a = l, b = b and b = 0. (4.8) 

(There is a second way of satisfying the NBP for the Kantowski-Sachs metric 
[§], but it will not lead to a universe containing black holes.) The path along 
the r Re -axis describes half of the Euclidean Nariai instanton S 2 x S 2 . Both 2- 
spheres have the radius H" 1 . Continuing parallel to the r Im -axis, the solution 
remains real: 



a(r Im ) 



-Re_ _2_ 
2H 



H- 1 cosh Hr lm , 6(r Im ) 



~2H 



(4.9) 



This describes half of a Lorentzian Nariai universe. Its spacelike sections 
can be visualized as 3-spheres of radius a with a "hole" of radius b punched 
through the North and South pole. This gives them the topology of S 1 x S 2 . 
Their physical interpretation is that of 3-spheres containing two black holes 
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at opposite ends. The black holes have the radius b and accelerate away 
from each other as a grows. The Nariai universe is a degenerate case of 
the Schwarzschild-de Sitter spacetime, with the black hole horizon and the 
cosmological horizon having equal radius J7j. 

The above path corresponds to half of a 2-sphere joined to a two- 
dimensional hyperboloid at its minimum radius H -1 , cross a 2-sphere of con- 
stant radius H~ x . It can be matched to any ag M > but only to bgM = H^ 1 
so the wavefunction will be highly peaked around that value of b. 

The first term of equation (|4.2|) vanishes and so the real part of the action 
for the Nariai solution comes entirely from the second term: 

= ~< = -J- (4-10) 

Now we compare the probability measures corresponding to the de Sitter 
and Nariai solutions. We find that in these models with a fixed cosmological 
constant primordial black holes are strongly suppressed, unless A is at least 
of order 1 in Planck units: 

exp (-2/* a c riai ) = exp (^j « exp = exp (-2/ d R c c Sittcr ) . (4.11) 

5 An Inflationary Model Without Black Holes 

Of course, we know that A 0, and therefore the models of the previous sec- 
tion are rather unrealistic. However, in inflationary cosmology it is assumed 
that the very early universe underwent a period of exponential expansion. It 
has proven very successful to model this behaviour by introducing a massive 
scalar field $ with a potential |m 2( l> 2 . If this field is sufficiently far from 
equilibrium at the beginning of the universe, the corresponding energy den- 
sity acts like a cosmological constant until the field has reached its minimum 
and starts oscillating. During this time the universe behaves much like the 
Lorentzian de Sitter or Nariai universes described above. 

But there are two important differences due to the time dependence of the 
effective cosmological constant A e g : Firstly, for the solutions of the Euclidean 
Einstein equations in the complex r-plane one can no longer find a path 
on which the minisuperspace variables are always real. However, we shall 
see that it is possible to satisfy conditions i to Hi of Section ^| by choosing 
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appropriate complex initial values. Secondly, it will be found in the next 
section that the black hole radius b is no longer constant during inflation. 

In this section, we introduce the massive scalar field for the model cor- 
responding to de Sitter spacetime, where the spacelike slices are 3-spheres 
containing no black holes. This model was first put forward by Hawking [§]. 
From the fluctuations in the cosmic microwave background as measured by 
COBE |J it follows that m is small compared to the Planck mass [fLOfl : 

m^0(l(r 5 ). (5.1) 

We will find complex solutions and the complex initial value of the scalar 
field, and we calculate the real part of the action. This has been done before 
by Lyons [|TTJ, but his paper contains a logical error to which we will come 
back later. 

The ansatz for the Euclidean metric is again 

ds 2 = N(r) 2 dT 2 + a(r) 2 dnl (5.2) 

Using the rescaled field 

cj) 2 = 4tt$ 2 (5.3) 

we obtain the Euclidean action 

377 r Ndra ( ^- + 1 - -a 2 ^- - -a 2 m 2 A , (5.4) 



4 J \N 2 3 N 2 

so that the effective cosmological constant is 

A eff (r) = m 2 0(r) 2 . (5.5) 
In analogy to equation (|3.3| ) we define 



H(r) = ^Ml) = (5.6) 

Variation with respect to a, <p arid iV gives the Euclidean equations of motion 
and the Hamiltonian constraint: 

- + l<P 2 + \m 2 <p 2 = (5.7) 
a 3 3 

+ 3-0-m 2 = (5.8) 
a 

- 2 -\-U 2 + \m 2 ct> 2 = 0. (5.9) 
a z a z 6 3 
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To evaluate $o(agM, 4>bm) using a semiclassical approximation we must find 
solutions in the complex r-plane that meet conditions i to in of Section |[ 
In particular, the NBP must be satisfied: 

a = 0, a = l, = O an d = for r = 0. (5.10) 

Assume that the initial value of the scalar field is large and nearly real: 

0*° » 1 » 00™ (5.11) 

An approximate solution near the origin is given by 

, , 1 



H Re^ H e r 

oo -r 
00 +E -,lnT n 



n=l 



(5.12) 
(5.13) 



for \r\ < O (l/# Re ) 



where the Taylor series is obtained by solving equation ( |5.8|) iteratively for 
0, using the NBP conditions ( |5.10| ) and the approximation ( |5.12| ) for a. It 
has the property that 

72«.+i = for all n. (5.14) 

We call equations (|5.12|) and (|5.13|) the inner approximation. Writing down 
the Taylor expansion explicitly to lowest non-trivial order 



l + A(#or) 2 O.V 
80o 



shows that is almost constant near the origin. 
As an outer approximation we use: 



(5.15) 



ao{r) 



im 



ipo + —^r + xo exp(3iH r) 



im 



a exp 
for < r Im < 



{T')dr' 



+ c exp 



im 



m 



(r')rfr' 



(5.16) 
(5.17) 
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While this solution does not satisfy the NBP, it will be good outside the 
validity of the inner approximation. Both the Xo-term and the co-term can 
be neglected for r Im ^> l/if Re , but they are useful for matching a and 4>o 
to a% and 0j at some |r| « 0(1/ H^ e ). Comparison with equation ( |5.15| ) 



shows that 



Ao 



O 



Im(xo) ~ 0. 



(5.18) 



In the region of the inner approximation, a will be nearly real on the 



Lorentzian line r Re 



2H* C 



Matching ao to aj fixes 



a 



2# Re 



c 



— i 
2F Rc 



(5.19) 



and ensures that a will remain nearly real on this line. To make </>(r) roughly 
real on the same line, by equations ( |5.16| ) and ( |5.18| ) we have to choose 



Im 



7T 



Re 



(5.20) 



in the outer approximation. 

o m in turn is fixed by matching <px to (f>o- Since it is very small, this 
requires evaluation of equation ( |5.13| ) to a very high order n. However, we 
need not calculate any coefficients since, by equation ( |5.14|) , ^Jf 1 is constant 
along the imaginary axis to any order n: 



Jm/ Im 



(5.21) 



Therefore it is convenient to choose a matching point tm on the imaginary 
axis: 



Re 



0, r 



Im 

M 



O (l/tf Re ) 



By equations ( |5. 16| ) and (|5.18| ) 4> 1 q is also constant along this axis: 



tlm/'^Im 
'O 



(r im ) 



Im 



(5.22) 



(5.23) 



so the result of the matching analysis will be independent of the precise choice 
of tm on the it should be. The matching condition is 



<P 1 o(tm) 



(5.24) 
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and by equations ( |5.20| ), ( |5.21| ) and ( |5.23| ) we obtain 



* = «° = ~w (5 ' 25) 

This result is non-trivial (e. g. ^ 'ipo' 6 )- We now see why the correct value 
for 0Q m is obtained in ref. |nj, although actually only ip 1 ™ is calculated there. 

We have thus satisfied condition ii of Section |2|. By the continuity of 
the outer approximation, condition in can be satisfied by fine-tuning </>Q m . 
Condition i is satisfied by the construction of the inner approximation. The 
only freedom left is the choice of 0^ e . This variable parametrizes the set of 
solutions. 

To calculate the Euclidean action for the solutions given above, we con- 
sider a path going along the real r-axis from the origin to r Re = and 
then parallel to the imaginary r-axis to tqm- Both a and are nearly real 
on the Lorentzian segment of this path, so the real part of the action can 
be approximated by an integral only over the first segment, using the inner 
approximation [|TT|]: 

# *4/; /2 " f l) « (5.26) 

The outer approximation is not valid after inflation ends, when 0^0. 
However, at this point we are already well inside the classical regime. A dust 
phase will ensue where oscillates; a and <fi will both remain real. Approxi- 
mate solutions for this regime have been given by Hawking and Page ||12|| . 



6 An Inflationary Model With Black Holes 

We now introduce a massive scalar field on a universe with spacelike sections 
S 1 x S 2 . Thus we will obtain a cosmological model similar to the Nariai 
universe of Section We find the complex solutions, initial conditions and 
the action in analogy to the previous section, but point out a few differences. 
Again we use the Kantowski-Sachs metric 

ds 2 = N(r) 2 dT 2 + a(r) 2 dx 2 + b(r) 2 dn 2 2 (6.1) 

and the rescaled field 

(j) 2 = 47r$ 2 . (6.2) 
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The Euclidean action is 



i «r , / b 2 b ab , 9 d> 2 ,., . 



2_2 i2 



+ 7T 



-d& 2 — 2a55 



T=0 



(6.3) 



and like in the previous section the effective cosmological constant is given 
by 

A eff (r) = m 2 0(r) 2 . (6.4) 
In analogy to equation ([O]) we define 



H(t)= ^A cfr (r) =m0(r). 



(6.5) 



Variation with respect to a, 6, and N gives the Euclidean equations of 
motion and the Hamiltonian constraint: 



b ab ^ ^2 

b ab 

b ab a - 9 9 , 9 
7 + ^ + -+0 2 + m 2 2 
b ab a 

0+ ( - + 2^ ) 0-m 2 
\ a o / 









(6.6) 
(6.7) 
(6.8) 



ab 



1 



2 1 

ab b 2 b 2 



\ 2 + mV 



0. (6.9) 
The NBP conditions corresponding to an instanton of topology S 2 x S 2 are: 



a = 0, a = 1, 6 = 6o, 6 = 0, = O and = for r = 0. (6.10) 

With the new definition (|6.5|) of if the inner approximation is given by: 

1 



ax(r) 
Mr) 



H Re^ H o e r 

OO 1 

+ E -f7nT* 



n=l 



nl 



m(px(T 



(6.11) 
(6.12) 
(6.13) 



for |r| < 0(1/^0 



rRe 
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The outer approximation is: 



<PO\J) 

a {r) 
bo(r) 



•00 + imf + Xo exp(iH r) 



ao exp —im 
1 

TO0 o (r) 
for < r Im 



(rOrfr' 



c exp 



im I (p(T')dr' 



(6.14) 
(6.15) 

(6.16) 



A matching analysis completely analogous to that of the previous section 
shows that a, b and will be nearly real on the Lorentzian line r Re 
if we choose the following initial values: 



71 



20, 



Re • 



m0 o 



(6.17) 



^ e is a free parameter. 



An interesting feature of the outer approximation is that the black hole 
radius grows with the horizon scale during inflation. On the Lorentzian line 
T Rc = ^h^- the field decreases linearly with time until it reaches zero and 
inflation ends. By equations (|6.14| ) and (|6.16|) b becomes very large on the 
timescale 



At, 



growth 



m 



(6.18) 



Again the inner approximation is used to calculate the real part of the 
Euclidean action. As in Section f| it comes entirely from the r = term: 



Re 



S 1 xS 2 



-7T 



bU e )' 



7T 



m 2 (0o e ) 



(6.19) 



7 The Probability for Primordial Black Holes 

In the previous two sections we have calculated the action for two inflationary 
universes. We now compare the corresponding probability measures 

= exp (^Wf) and Ps ' M ^ ] = exp fcw ) • (7 - 1) 
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The universe containing black holes is heavily suppressed, if 0q * s n °t large 
enough to make the initial effective cosmological constant equal to the Planck 
value. Thus the formation of black holes with initial sizes significantly larger 
than the Planck scale is very unlikely. The semi-classical approximation 
should be good in these situations, so one can have confidence in this con- 
clusion. 

The semi-classical approximation will break down for solutions with initial 
cosmological constants of the Planck value in a region where the curvature is 
on the Planck scale. However, this region contributes an action less than one 
in Planck units and one would not expect quantum effects to change this. 
Thus it seems clear that the only primordial black holes with any significant 
probability start with no more than the Planck size: 



The Nariai solution is unstable to quantum fluctuations [7|]. At the begin- 
ning of inflation it becomes a non-degenerate Schwarzschild-de Sitter space- 
time. Once the black hole horizon is inside the cosmological horizon the 
black hole will start to lose mass due to Hawking radiation. If the black 
hole horizon is somewhat smaller than the cosmological horizon, the black 
hole will evaporate and disappear. However, there is a significant probability 
that the areas of the two horizons will be nearly enough equal for them to 
increase together. The consequences of this result for the global structure of 
the universe will be presented in a forthcoming paper. 
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